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Abstract. 

We study the non-local superconducting pairing of two interacting Anderson impurities, 
which has an instability near the quantum critical point from the competition between the 
Kondo effect and an antiferromagnetic inter-impurity spin exchange interaction. As revealed 
by the dynamics over the whole energy range, the superconducting pairing fluctuations acquire 
considerable strength from an energy scale much higher than the characteristic spin fluctuation 
scale while the low energy behaviors follow those of the staggered spin susceptibility. We argue 
that the glue to the superconducting pairing is not the spin fluctuations, but rather the effective 
Coulomb interaction. On the other hand, critical spin fluctuations in the vicinity of quantum 
criticality are also crucial to a superconducting pairing instability, by preventing a Fermi liquid 
fixed point being reached to keep the superconducting pairing fluctuations finite at low energies. 
A superconducting order, to reduce the accumulated entropy carried by the critical degrees of 
freedom, may arise favorably from this instability. 



1. Introduction 

In conventional superconductors, the electron pairing arises from the mediation of lattice 
vibrations [lj. In cuprate superconductors, the pairing mechanism for superconductivity has 
been hotly debated as to whether the pairing occurs via the coupling between fermionic 
quasiparticles and bosonic modes, or whether a bosonic glue is really necessary [2]. In the 
former type of mechanism, the scenario of electronic coupling to spin fluctuations is the leading 
contender [3j. However, in high-T c cuprates, the antiferromagnetic Neel state is terminated 
far before the unconventional superconductivity emerges as the hole doping varies. As such, 
the relation between a magnetic quantum criticality point and high-T c superconductivity has 
not been developed. Furthermore, the final identification of the mechanism is also complicated 
by the non-negligible electron-phonon coupling as revealed in various experiments [I]. In this 
aspect, heavy fermion materials have been providing test grounds for novel mechanism for 
superconductivity [5]. In these systems, the interplay between magnetism and superconductivity 
has been observed unambiguously. The generic phase diagram shows that, when a tuning 
parameter such as doping, magnetic field, and pressure is varied, the magnetic state is either 
interrupted by or in close vicinity to the appearance of a superconducting state. Therefore, they 
are the most ideal systems for the study on unconventional superconductivity associated with 
strong electron correlations in general and the quantum criticality in specific. 



To address the relation between superconductivity and magnetic quantum criticality, we 
consider here a two-impurity Anderson model, which is a minimal setting for antiferromagnetism 
and non-local superconducting pairing. It is also the quantum impurity/cluster model for the 
periodical Anderson lattice model or the Hubbard model within the cluster dynamical mean-field 
theoretical approach [6]. This model captures the essence of the competition between the local 
Kondo coupling between spins of the localized electrons to a conduction band, and the non-local 
spin exchange interaction between spins of localized electrons, e.g., Ruderman-Kittel-Kasuya- 
Yosida (RKKY) interaction in heavy fermion materials. Due to this competition, a quantum 
phase transition occurs between a Kondo resonance state and an inter-impurity spin singlet 
state [3 El O [TU1 HH Q21 Q21 E], and multiple energy scales can show up [TH]. In this paper, we 
discuss the superconducting pairing between two Anderson orbitals and its mechanism. 

2. Model 

The Hamiltonian for the two-impurity Anderson model can be written as 

H = E ^cL^ + E( Fkeikr<c L/- + hx -) + E( e //*U- + f + /s i • s 2 >(!) 

k<r kcri ia 

which describes two interacting local orbitals fi a (i = 1,2), with energy level ej and onsite 
Coulomb repulsion U, in hybridization with a non- interacting conduction electron band 
with the strength Vk at each impurity site rj. In terms of the even or odd parity combinations 
fp=(e,o),(7 = (/lo- =t /2<t)A/2, the fluctuations due to conduction electrons can be represented by 
two separate baths with the hybridization functions r ej0 (u;) = — ^Im[^ k V^|e jk ' ri ±e* k ' r2 | 2 / (uj — 
e^ + iO + )]. In a realistic system RKKY interaction is perturbatively generated. Here we consider 
a case T efi (uj) = Tq for \uj\ < D (To is a constant and D the half bandwidth of the conduction 
electrons), for which the generated RKKY interaction vanishes, and we tune the intersite spin 
interaction by a direct RKKY term /Si • S2. 

We focus on a non-local superconducting pairing correlation function Xsc = ((A sA)) with 
At = (fl^-fl^ + f^f\_^)/\/2, i- e -> ™ a singlet configuration. To illustrate its relation to spin and 
charge fluctuations, we also provide results on the uniform and staggered charge susceptibilities 
Xc(u,a) = (( n i =t n2',ni ± na))/2, the uniform and staggered spin susceptibilities \u,a = 
((S lz ±S2 z ;S lz ±S 2z )}/2, the current fluctuation Xcurr with k\ urr = iY^aiflaha- ft a ha)/V^, 
as well as the single-particle properties including the spectral function Af = — (l/7r)ImG/ 
and the interaction part of the self energy — ImSj. These quantities are calculated from the 
complete-Fock-space numerical renormalization group method (CFS-NRG) [16J while a two- 
particle Green's function method [17] is adopted for calculating the self energy. 

3. Correlation fluctuations and the associated energy scales 

The results for these dynamical quantities for different values of I are shown in Fig. [TJ Some 
of the calculation details can be found in Ref. [15]. With the tuning of an antiferromagnetic 
/, there is a change of the ground state from the Kondo resonance state Aj(uj = 0) « l/(7iTo) 
[e.g., Fig. QJa-b)], to the inter-impurity spin singlet state with a pseudogap Af(u>) ~ uj 2 [e.g., 
Fig. HJd)]. Though the spectral weight has a sudden change at the Fermi energy, this transition 
is continuous as the effective (local) Fermi liquid (FL) temperature Tp (both competing states 
belong to FL fixed points) vanishes uniformly at the critical point I c ~ 2.3T/< [Fig. [TJc)], where 
Tk is the single-impurity Kondo temperature. \ u and Xa have different behaviors. While x"( w ) 
can be characterized by one monotonically increasing energy scale T s f, Xa( w ) m general has 
a flattened peak spanning a finite energy range from Tp to T s f. (Tp and T s f correspond to 
Tl and T# in Ref. |15|). T s f therefore characterizes a spin fluctuation scale, which takes the 
value between Tk and I. Between Tp and T s f, there are still finite spectral weights but these 




Figure 1. Dynamical properties of the two-impurity Anderson model for various values of the 
inter- impurity spin exchange interaction I: (a) 1 = 0, (b) I < I c , (c) I ~ I c and (d) I > I c . 
The spectral functions Aj (in unit of D), and the imaginary part of the self-energy — ImEj 
(in unit of 1/D) are shown in log- linear scale (left y-axis), while the imaginary parts of the 
uniform spin susceptibility x'ui the staggered spin susceptibility x'Li the uniform and staggered 
charge susceptibilities (the same in this case) Xc> an d the superconducting pairing fluctuation 
Xsc (also the half of the current fluctuation Xcurr due to the axial charge symmetry), all in unit 
of 1/Tk, are shown in log-log scale (right y-axis). Different energy scales, such as the Fermi 
liquid scale Tp, the spin-fluctuation scale T s f, and the local moment scale Tpm are labeled. Here 
Te.o = 0.045-7T.D and ej = —U/2 = —D ~ — 10 3 Tx- In this case, Af and — ImSj are found to be 
the same for different parities and spins. 



quasiparticles develop non- Fermi liquid behaviors, which relate to partially screened Kondo 
resonances |15j . 

We now examine other correlation functions. For the charge fluctuations, we find that they 
play little role in the low energy physics. They have resonance peaks around U/2 and decrease as 
~ oj at low energies from a scale Tlm- Here, the uniform and staggered charge susceptibilities are 
found to be the same and remain unchanged with the tuning of /. (However, in the particle-hole 
asymmetric case, it is found that the uniform charge susceptibility also diverges at the quantum 
critical point [18j QjJ]). For the superconducting (SC) pairing correlation function Xsc(^>), we 
find the low energy part of Xsd^) follows the behavior of Xa( u; )> as ~ w for cj < Tp and as ~ C sc 
for oj > Tp, where C sc is relatively a constant. At the critical point, as Tp — > 0, x'sci®) becomes 
divergent, leading to a SC pairing instability. However, the high energy part doesn't follow the 
spin fluctuations: it has a resonance peak at U and the SC fluctuation softening happens from 



Tlm, which is much larger than T s f, down to Tp, or Xs C (0) ~ C sc ]xl(Tlm /Tf)- 

Tlm is the also the energy where — ImS(w) reaches a maximal value. In the single-impurity 
Anderson model [20J, the renormalization group (RG) flow starts from the free orbital fixed 
point, rather than directly from the local moment fixed point as in the Kondo model. The free 
orbital fixed point is unstable against any interaction U, or the hybridization with the conduction 
electrons T. U e ff and T e ff increase at first, approaching the local moment fixed point, which 
is characterized by U e ff — > oo and a small Kondo coupling Jic,eff = ^eff/U e ff, and then U e ff 
decreases while r e ff keeps increasing, reaching the strong coupling fixed point Jx,eff oo (due 
to the marginal nature of Jk, such change is slow and only becomes abrupt near Tlm 
therefore characterizes the proximity to the local moment fixed point in the RG flow. 

4. On the superconducting pairing 

The fact that the SC pairing fluctuations become significant before either Tk or /, and in 
addition their weight remains relatively unchanged at uj = T s f for finite 7s, implies that the 
SC pairing might not require the spin fluctuations as the glue. The central issue is therefore to 
understand the SC pairing fluctuation softening beginning at Tlm- I n terms of the even and 
odd parity orbitals, the interaction part of the original Hamiltonian can be rewritten as 

Hf, int = (17/2 - 37 /8)(n et n el + n ot n ol ) + (U/2 + 37 /8)J e • J - (U/2- I /8)S e • S G , (2) 

where J and S are axial charge [7] and spin of the even and odd orbitals. For J + = fl^-fpp it 
can be easily shown that the transverse "staggered" and "uniform" axial charge susceptibilities, 
((J£ "F J^'i^e "F Jo)) correspond to the non-local singlet SC pairing correlation function 
and the local SC pairing correlation function xfoc = ((Hi fi-tfJi'iHi filfit)) ■ Due to the 
axial charge rotational symmetry, they are the same as twice of the correlation functions of 
Jg =p = (n e =Fn )/2, which relate to the current fluctuations and uniform charge susceptibility, 
respectively (explicitly verified by calculations). With an "antiferromagnetic" axial charge 
exchange interaction I ac = U + 37/4, the staggered component has negative energy while the 
uniform component has positive energy, i.e, the non-local SC pairing is energetically favored 
against the local SC pairing. Such axial charge fluctuations are also subjected to the coupling 
to the conduction electron bath. In Fig. [2j we show the behaviors of x'sc an d (xfoc)" ^ or different 
values of the hybridization constant Tq (U is fixed). We observe that the non-local SC pairing 
fluctuations are always stronger than the local ones, indicating the "antiferromagnetic" nature. 
Besides, the enhancement (or the softening) of x'sc 1S only achieved by a relative large value 
of Tq. While a detailed analysis of this behavior is desired, we can associate it with an "axial 
charge Kondo" effect, by making an analog between the axial charge and the spin. We expect 
a "Kondo" -like exchange interaction between the axial charges of the local orbitals and those of 
the conduction electrons in the form J ac J p - J c ,p, and as expected, J ac ~ Jk- Such an interaction 
competes with I ac in the same fashion as the competition between the Kondo coupling and 
RKKY interaction in the spin channel. Indeed, we find the behaviors of x'sc m Fig- [2] share 
similarities to the spin susceptibilities when the ratio between Tk and 7 is varied. In the spin 
dynamics, the softening of x'L ls due to the equivalent strength of local Kondo splin-flip scattering 
and inter-impurity spin singlet to triplet scattering. The same type of physics is expected for 
the axial charge fluctuations: the axial charge fluctuations with conduction electrons compete 
with the axial charge fluctuations between non-local ("singlet") and local ("triplet") SC pairs. 

On the other hand, the low energy edge of the SC pairing fluctuation softening is Tp, rather 
than T s j for finite 7s. This indicates that the SC pairing fluctuations are only suppressed when 
the strong coupling fixed point, or a FL fixed point is reached. This can be understood as that the 
effective Kondo couplings in the spin and axial charge channels are the same. When the strong 
coupling fixed point is reached, Jx,eff ~^ o°; the axial charge for each parity forms an "axial 
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Figure 2. The non-local singlet 
iXsc-, solid lines) and the local (xfo C > 
dashed lines) superconducting pairing 
fluctuations for different values of 
T efi (uj) = Tq. U is fixed as 2D, as 
in Fig. 1, and 1 = 0. 



charge Kondo singlet" with those of the conduction electrons. The axial charge fluctuations 
between even and odd parities, i.e., the SC pairing fluctuations are suppressed. While the spin 
exchange interaction / (antiferromagnetic in our study, but we expect similar physics for the 
triplet pairing with ferromagnetic /) can create an "antiferromagnetic" interaction (I a c) m the 
axial charge channel as well [Cf. Eq.([2])], but in general, I <^U. This also explains that x'Li 1 ^) 
has little change at u = I. The effect of / in this picture is rather to compete with the Kondo 
effect and keep Jx,eff finite at low energies. This implies that the SC pairing can only keep its 
strength to low energies when the Fermi liquid fixed point is avoided, i.e., near the quantum 
critical point. This can also be evidenced from a low energy effective Hamiltonian in terms of 
Fermi liquid parameters [7j, where I a c,eff is found to be divergent at the quantum critical point. 

In the traditional theory of superconductivity, the formation of Cooper pairs requires an 
attractive interaction, direct or generated with certain glues. Our results rather imply that such 
an effective attractive interaction in the non-local singlet channel is originated from the onsite 
repulsive interaction. The SC pairing fluctuations at high energies may simply be two-particle 
excitation fluctuations. The "antiferromagnetic" nature of I ac , which appears as an attractive 
interaction in the site basis — (I ac/ 2) ( n it n 2i + n 2t n u)> m &y be due to the repulsive nature of 
local Coulomb interaction, i.e., the non-local two-particle excitations are favored against local 
pairs. The central issue revealed by this study is that such two-particle excitations can be 
preserved to low energies. This is achieved by an "axial charge Kondo effect" to enhance pair 
fluctuations with coupling to conduction electrons. Such an enhancement maintains further 
to the Fermi energy when the spin dynamics are critical near a quantum critical point. We 
notice that the spin fluctuations are the dominating fluctuations at low energies. In other 
words, the weights of the low energy states are dominated by the single-occupancy states (in 
the site basis). The SC pairing fluctuations are subjected to the weight of the double or zero 
occupancy configurations in the low energy states, which is ~ (T/U) 2 . Therefore, the effective 
SC pairing energy scale is I a c,eff = U(T/U) 2 = T 2 /U instead of U . The dominating excitations 
would still be single-particle excitations. However, the single-particle excitations at quantum 
criticality have non-Fermi liquid behaviors. In the two-impurity model, they are formulated as 
Majorana Fermions |1Q[ 111} [T3]. The origin for Majorana Fermions is also related to the finite 
JK,effi which can only lead to partial-screened Kondo resonances. They carry finite residue 
entropies [10j [Til 03] • A superconducting order, by forming Cooper pairs which is a lower 
entropy state, becomes favorable to these single-particle excitations. We notice that the entropy 
accumulation is also a generic thermodynamic feature of quantum criticality |21| . 

In summary, from the dynamical property of the superconducting pairing correlation function 
in the two-impurity Anderson model, we learn that the superconducting pairing doesn't 



necessarily need spin fluctuations, or bosonic modes as the glue. It can be originated from 
the effective Coulomb interaction. However, such pairing strength at low energies only remains 
finite when spin fluctuations are critical to avoid a Fermi liquid fixed point. This is manifested 
as being in the vicinity of certain quantum critical points. This provides us a new picture on 
the origin of unconventional superconductivity in strongly correlated electron systems. It will 
be interesting to check how this idea can be applied to the lattice systems. 
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